Dynamics of formation and decay of coherence in a polariton condensate 
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We study the dynamics of formation and decay of a condensate of microcavity polaritons. We 
investigate the relationship between the number of particles, the emission linewidth and its degree 
of linear polarization which serves as the order parameter. Tracking the condensate formation, we 
show that coherence is not determined only by occupation of the ground state, bringing new insights 
into the determining factors for Bose-Einstein condensation. 



Whereas a lot is known of the thermodynamics of 
quantum fluids — at equilibrium and in infinite size 
systems — the dynamics of formation of a Bose-Einstein 
condensate (BEC) is still vastly an open question. This 
dynamics is not easily detectable in atomic BEC and 
has not been studied experimentally, due to the very 
short timescale in which cold atomic gases reach ther- 
mal equilibrium. A good alternative to investigate the 
process of condensate formation is that of microcavity- 
polaritons [H, although this system differs from atomic 
condensates due to its intrinsically out-of-equilibrium 
character. Polaritons arise from the strong coupling of 
photons and electrons in semiconductor microcavities Q 
and have recently shown to undergo a non-equilibrium 
phase transition with a spontaneous buildup of coherence 



in the ground state [1,0,1 
observed in atomic BEC 



similar to what has been 
|. Thanks to the small po- 
lariton lifetime, various properties can be tracked contin- 
uously from the photons they emit. This has been used 
to measure the spatial first order coherence jl0[] , the tem- 
poral first and second order correlation functions [ill. 12| 
and even real time dynamics such as propagation [131. 

The dynamical evolution and origin of the degree of 
coherence in a polariton condensate remain however un- 
clear. Intriguing results have been recently reported of 
one order of magnitude increase of the ground state po- 
laritons coherence time, by reducing particle fluctuations 
in the excited states [12]. Such experiments reveal the 
possibilities for the dynamical observation of condensa- 
tion [14]. Also some theoretical studies have been re- 
ported based on stochastic simulations of the order pa- 
rameter jTjjj. In this work, we pursue this goal both ex- 
perimentally and theoretically, and track the dynamics 
of formation of a polariton BEC, following how the con- 
densation forms and decays in a pulsed experiment. 

Since polaritons have a spin degree of freedom that 
is passed to the polarization of the photon they emit, 
the buildup of off-diagonal elements in the density ma- 
trix (characteristic of coherent states) can be accessed 



through the degree of linear polarization: 

A = 2 | S^oora^) |/«4 T aoT> + <4| tt o|» , (1) 

with aoT and aoj, the field operators for spin-up (right- 
circular polarization) and spin-down (left-circular polar- 
ization) polaritons, respectively. If spin- up and spin- 
down polaritons are uncorrelated, with a thermal distri- 
bution and no determined phase relationship, the light 
is unpolarized. If, however, the two spin-projections are 
condensed with well denned order parameters (aoj), (a>oi) 
that factorize the numerator of Eq. (pQ), the emitted light 
acquires a definite polarization. Therefore, D\ is linked 
to the degree of coherence of the condensate and to the 
imbalance of spin- up/spin-down polaritons [6l.Hl. lrj 17j |. 
The spin degeneracy can be lifted producing a Zeeman 
splitting of polariton energies. If we can further assume 
that the two spin populations have similar evolutions, 
these can be modelled with a convolution of thermal 
and coherent states 1(3, 18[ with time-dependent thermal 
and coherent fractions. In our case, we have observed 
experimentally that, except for a short transient after 
the pulsed excitation, spin-up and spin-down populations 
equalize and evolve similarly in time Therefore the 
degree of linear polarization is directly related to the de- 
gree of second order coherence (at zero delay) of the 
polariton gas: D\ 



7 (2) 



If the spin populations 



have very different dynamics, one has to keep track of the 
coherence degree of each of the two spin components and 
this simple expression no longer holds. Up to now, this 
relationship has been exploited for measurin g th e degree 
of coherence only in stationary situations [6Ull| . In this 
work, we use it to study experimentally the dynamics of 
condensation. Using simultaneous detection of energy- 
and time-resolved photoluminescence, together with the 
degree of linear polarization of the ground state, we mon- 
itor the formation and decay of the polariton condensed 
phase in the ground state after a circularly-polarized, 
pulsed excitation at high energy. Comparing our data 
with a phenomenological model, we investigate the evo- 



2 



Energy (eV) 



o 

50 
100 
3150 
E 200 
250 
300 



I 




Energy (eV) 




200 300 
Time (ps) 

FIG. 1: (Color online) Experimental results on the dynam- 
ics of BEC formation after the arrival of the non-resonant 
pulse, (a) Emission energy from the ground state (k = 0) as 
a function of time (left) and the corresponding photolumines- 
cence intensity (right). The color code goes from black, no 
emission, to red, intense emission, saturating in a black cen- 
tral region, (b) Spectral lineshapes extracted from (a), (c) 
In blue dots, average condensate population no (normalized 
photoluminescence intensity); in purple squares, linewidth ko 
(meV) of the emission peak, limited by the spectral resolution 
of ~ 0.1 meV; in brown rhombus, degree of linear polarization 
D\. Times for coherence build-up (£i), maximum population 
(£2) and coherence (£3) are marked with dashed lines. Inset: 
polarization polar plot at £ = £3 with maximum coherence for 
the present case (I) and for a second point in the sample (II). 



lution of the condensate statistics and coherence. 

In the experiments, polaritons are created in a CdTe- 
based microcavity grown by molecular beam epitaxy. 
The sample consists of a Cdo.4Mgo.6Te 2A-cavity with 
top (bottom) distributed Bragg reflectors of 17.5 (23) 
pairs of alternating A/4-thick layers of Cdo.4Mgo.6Te and 
Cdo.75Mno.25Te. Four CdTe QWs of 50 A thickness, sep- 
arated by a 60 A barrier of Cdo.4Mgo.6Te, are placed 
in each of the antinodes of the electromagnetic field. 
Exciton-photon coupling is achieved with a Rabi split- 
ting of 23 meV; the sample is kept at 4 K and condensa- 
tion is obtained around zero detuning between the exci- 
ton and the bare cavity dispersion. Photons are injected 
via circularly polarized, ps-long pulses tuned at the first 
minimum above the Bragg mirror stopband (~ 40 meV 
from the bottom of the dispersion) in order to guarantee 
the memory loss of the laser polarization and the coher- 
ence by the time polaritons are formed at the bottom 



of the band. Light escaping from the microcavity is dis- 
persed through a 0.5 m spectrometer and analyzed with 
a streak camera. This allows to resolve in time the in- 
tensity and energy of the emission of the ground state 
(FigOJa) and (b)), and simultaneously analyze the direc- 
tion and linear degree of polarization of the emitted light 
together with the polariton linewidth (FigQJc) and in- 
sets). The direction of the linear polarization is pinned to 
a crystallographic axis, as has been observed with contin- 
uous pumping in a sample of the same batch growth 111! 
and with pulsed pumping in samples similar to ours [20[ . 
Without this pinning, the polarization would average to 
zero, which is detrimental to our approach but is other- 
wise one of the most direct demonstration of symmetry 
breaking associated to the quantum phase transition, as 
reported by Baumberg et al [6]. In our case, the direc- 
tion of the pinning depends on the point of the sample, 
as shown in the insets (I) and (II) of FigQJc). 

The non-resonant laser excitation creates a large pop- 
ulation of carriers, at high energies, that form the polari- 
tons. Thanks to the quick relaxation mechanisms, like 
phonon emission and polariton-polariton scattering [211 ] . 
the population of polaritons relaxes into the fundamental 
state (with in-plane momentum k = 0). Figure [T] depicts 
the experimental time evolution of the ground state after 
the arrival of the pulse at some time £ < 0, till the com- 
plete decay of the condensate. The intensity of the light 
emitted by the condensate, first shown in Fig. []Ja)-right 
(in density plot), is normalized to represent the ground 
state population no in Fig. [2(c) (blue circles). The abrupt 
change into a nonlinear growth — associated with the on- 
set of stimulated enhancement — is matched to 1. We 
note from the ground state energy (Fig. Ufa)), that the 
dispersion is already blueshifted (as compared to the long 
time energy, at £ > 300 ps), even when no is very small 
(at £ ~ 0). This confirms that the blueshift is mainly 
due to exciton screenin g an d carrier-carrier interaction 
in the exciton reservoir 22|, 23, 0], and that it is only 
slightly affected by the ground state occupation. Note 
that the blueshift is always much smaller (< 1 meV) 
than the Rabi splitting, so that the system remains in 
strong-coupling. As the reservoir depletes with time, the 
overall population and the blueshift monotonically de- 
crease, while no first increases and then decreases. The 
coherence imprinted by the laser pulse and its original po- 
larization is lost during this efficient cascade of decays, 
populating the ground state with initially uncorrelated 
polaritons. This is clear from the linear polarization D\ 
(brown rhombus in Fig. [2(c)), that takes some time to 
build up after no ~ 1. Passed this particular time (corre- 
sponding to £1 = 61 ps in the figures) the system is above 
threshold with a condensed fraction and a phase locking 
of the spin- up and spin-down condensates. Stimulated 
scattering into the ground state becomes the dominating 
process and the population no grows up to a maximum 
(with max(no) ~ 25 at £2 = 118 ps). However, the co- 
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FIG. 2: (Color online) Hysteresis in D\ (a) and kq (c) as a 
function of no, with time running in the sense of the arrows, 
(a) Three experiments are shown with increasing excitation 
powers (the one under study in brown rhombus) . (b) Growth 
with power of the delay between maximum population and co- 
herence for these three cases. The linewidth of Ref. [3], under 
CW excitation, is given in (c) (dashed green) for comparison. 



herence reaches its maximum, max(Di) = 0.75, at a later 
time, ts = 140 ps, with a proper dynamics that does not 
follow instantaneously that of no: D\ decays for another 
150 ps at a different rate than that of the population. 
D\ is represented as a function of no in Fig. [21(a) > and 
compared with two other cases in order to make evident 
the existence of an hysteresis loop. As the excitation 
power is decreased, not only no and D\ decrease, but the 
hysteresis also gradually disappears. 

The linewidth of the ground state emission is also dra- 
matically affected by the condensation. In the simplest 
picture of a single mode with an effective pumping rate 
from all the other states Wq u1 = Sfc^o(l + n k) (nk 
the population of the fcth state) and a radiative decay 
rate 70 at k = 0, the linewidth can be approximated by 
= (Wo ut + 7o)/(l + no). This expression is enough 
to understand qualitatively the linewidth dynamics. As 
shown in Fig. QJc) (purple squares), k,o drops dramat- 
ically at the phase transition due to the exponential 
growth of no- Close to the maximum population, the 
linewidth reaches the limiting value of our spectral reso- 
lution (0.1 meV). Note that ^0 not only depends on no 
but is also affected (broadened) by the populations of the 
levels feeding the ground state, which are enhanced be- 
fore the transition takes place. The steady state value of 
tzo — » 0.9 meV (t > 300 ps) corresponds to the linewidth 
of the bare ground state polaritons, 70. This results in 
another hysteresis shown in Fig. [21(c). Although polari- 
tons, being interacting bosons, also suffer screening due 



to Coulomb interaction, this effect is only significant at 
extremely high polariton densities (25[. This regime is 
never reached in our experiment. Kasprzak et al. [3[ ob- 
served under CW excitation an increase in the linewidth 
(dashed, green line in Fig. [21(c)) and attributed it to 
polariton-polariton interaction. However, we find that 
the ground state population is too small to produce this 
effect (the particle number should be three orders of mag- 
nitude higher). The reason for such a significant increase 
of the linewidth is more likely to be found in the high 
number of reservoir particles, much higher than no in the 
case of CW experiments. With our time dependent ex- 
citation, excluding the first few tenths of ps, the amount 
of particles in the reservoir is significantly smaller than 
that at the bottom of the polariton dispersion, reducing 
drastically this extra decoherence effect. 

Our experiments on polariton coherence formation are 
supported by a simple two-level model (sketched in the 
inset of Fig. [21(a)), that links together the dynamics of the 
coherence degree with the population and the linear po- 
larization. We work under the assumption — supported 
by the observation — that both spin components evolve 
similarly following this dynamics. We proceed in the 
spirit of the evaporative cooling description [26[ , that has 
been used before in models including polariton-polariton 
and photon-polariton scattering [HI, [I?!. A reservoir, 
state 1, is populated in an incoherent way from higher 
levels, what we describe with an effective time-dependent 
rate P\(t) which has a Gaussian profile in time, rem- 
iniscent of the excitation of the laser at very high en- 
ergy, with temporal effective width At pu i se . These par- 
ticles are scattered incoherently into the ground state, 
level 0, at a different rate Piq. Both levels loose po- 
laritons with decay rates 70 and 71, respectively. We 
also include the inverse process of promoting polaritons 
up the branch at some small rate 701. This dynamics 
is described with Lindblad terms in a master equation 
of the ground-state/reservoir system. This leads to rate 
equations for the time dependent distributions of parti- 
cles, which are solved numerically to compute the ob- 
servables of interest (no, g^ 2 \ D\, etc. . . ). The linewidth 
is obtained from ^0 ~ [701 (1 + ^1) + 7o]/(l + n o)- In 
this model, the coherence builds up spontaneously, while 
other works require a seed (such as an initial coher- 
ent state) in order to investigate the dynamics of coher- 
ence. The results, plotted in Fig. [U(a) with the same 
color code, are qualitatively similar to the experimental 
data in Fig. [U(c). With no need for more levels, we also 
reproduce the existence of a delay between maximum co- 
herence and population. We find that when the effective 
pulse P\(t) gets sharper in time (small 1/At pu i se ), not 
only max(no) and max(Z}i) are larger, but also the de- 
lay between them increases, as can be seen in Figs. [31(b) 
and (c). The delay can even become negative when the 
pulse is too flat. The origin and sign of the delay can be 
explained as follows: for a sharp and strong pulse, the 
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FIG. 3: (Color online) (a) Results from the theoretical model 
(sketched in inset): no (blue), D\ (brown) and (purple), as 
a function of time (rescaled to ps for comparison) . In (b) , the 
maximum max(no), max(Di) and minimum min(^o) values 
reached in the simulations as the effective pulse narrows (with 
increasing 1/A£ pu i se ) down to the case in (a). The delay £3 — £2 
between max (no) and max(Di) is also plotted in (c). 



of the dynamics of a polar iton condensate, created un- 
der a non-resonant, circularly polarized pulse, and the 
study of its coherence buildup and decay. We track the 
order parameter of the transition through the degree of 
linear polarization, that spontaneously builds up. This 
occurs when the average population of the ground state 
exceeds one and is accompanied by an abrupt decrease 
in the emission linewidth. The maximum degree of co- 
herence that is achieved does not coincide in time with 
the maximum population, evidencing that the coherence 
of the condensate has its own particular dynamics. The 
pulsed excitation regime is a worthy tool revealing neatly 
the interplay of all the relevant quantities of polariton 
condensates and displaying the mechanisms of formation 
and decay of a BEC. We support our claims with a simple 
theoretical model. 
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ground state population no rises so quickly that coher- 
ence formation cannot simultaneously build up (through 
stimulated interaction with the upper level), and settles 
to its own, intrinsic dynamics instead (the delay becomes 
independent of the pulse duration as this one gets nar- 
rower). With a flat pulse, on the other hand, no grows 
slowly and weakly, giving time to the condensate coher- 
ence to form together with it, possibly overtaking the 
dynamics for very long effective pulses and yielding a 
negative delay (Fig. E9(c)). Increasing the in- flow of par- 
ticles (e.g., increasing P10/701)? a l so nas a similar effect 
of enhancing the coherence and its rising time. On the 
other hand, the total intensity of the pulse is not deter- 
mining the delay. This argument can be extended to the 
experiment, where a more intense pulse also results in a 
sharper profile for the effective excitation, given that the 
chain of de-excitations towards level 1 is more efficient, 
thanks to the density dependence of the exciton-exciton 
scattering (2^ . 

Finally, although our simple model does not lead to a 
quantitative agreement with the experiments (time has 
been rescaled for a better comparison), it is of a funda- 
mental character, since it contains the minimum num- 
ber of elements to reproduce BEC formation and links 
together its essential features (population, coherence, 
linewidth and polarization). The model should include 
more levels in the de-excitation chain towards the ground 
state to reach to a quantitative agreement, not bringing, 
however, any essential new features into the description 
(as we have checked but which we do not discuss here). 

In conclusion, we report the experimental observation 
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